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Abstract 

We present the design of two-dimensional photonic crystal microcavities in 
thin diamond membranes well suited for coupling of color centers in diamond. 
By comparing simulated and ideal field distributions in Fourier and real space 
and by according modification of air hole positions and size, we optimize the 
cavity structure yielding high quality factors up to Q = 320000 with a modal 
volume of V e g = 0.35(A/n) 3 . Using the very same approach we also improve 
previous designs of a small modal volume microcavity in silicon, gaining a 
factor of 3 in cavity Q. In view of practical realization of photonic crystals 
in synthetic diamond films, it is necessary to investigate the influence of 
material absorption on the quality factor. We show that this influence can 
be predicted by a simple model, replacing time consuming simulations. 
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1. Introduction 

Microcavities in two-dimensional photonic crystal slabs allow to strongly 
confine light in volumes smaller than one cubic wavelength. They are ex- 
pected to enable the realization of single photon emitters PQ, low threshold 
nanolasers [2H3], ultra small filters [5] and highly-efficient emitters in which 
the spontaneous emission of light from single emitters is controlled at the 
quantum level [SHE]- Due to highly developed processing techniques for semi- 
conductor materials most of these quantum information devices have already 
been demonstrated in silicon or GaAs. Besides the use of semiconductor ma- 
terials for quantum information processing, diamond has attracted significant 
interest in recent years due to the extraordinary properties of optical active 
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defect centers that can be controlled at room-temperature [HI EI]- There have 
been a number of recent proposals for employment of color centers in dia- 
mond for cavity enhanced single photon sources [TT] , for cavity enhanced spin 
measurements [12] or as optical qubits in quantum networks [13], quantum 
gates [H] and measurement-based quantum computing [T51 HE] . 

All these proposals require direct coupling of an emitter to a cavity mode 
with high quality factor Q and small modal volume V e g. In such a cavity the 
spontaneous emission of an emitter, placed in the maximum of the electric 
field of the cavity mode, is enhanced by the Purcell- Factor F [TT] : 



In this weak coupling regime the cavity figure of merit scales as Q/V e e, 



whereas the dynamics of strong emitter-photon coupling scales as Q/yV e s 




In recent years, several photonic crystal cavity designs have been pro- 
posed, e.g. modified line- or point defects in two-dimensional photonic crystal 
membranes. In silicon, modified waveguides structures have been fabricated 
yielding very high quality factor of 10 5 - 2.5 ■ 10 6 [X9H2T] with modal volumes 
of 1.2 - 2(A/n) 3 . For diamond-based modified waveguide structures similar 
theoretical Q-factors with a modal volume of about V e s ~ 1.7(X/n) 3 have 
been predicted [22H21]- In practice, the cavity Q-factor in diamond might 
strongly be limited by optical loss in the diamond material [25] . 

In order to obtain large enhancement of the spontaneous emission rate 
for cavities with modest Q-factors, point defect cavity structures with even 
smaller modal volumes as compared to waveguide section cavities are re- 
quired. Different point defect structures with zero, one or more missing holes 
have been proposed [25]. As shown by Akahane et al. [27] the quality factor 
of a point defect cavity can be improved further by optimizing the surround- 
ing photonic crystal structure. They have demonstrated that the Q-factor 
of a photonic crystal in silicon with three missing holes (M3-cavity) can be 
increased up to 4.5 • 10 4 [27] by optimizing the next neighboring holes ac- 
cording to the method of "gentle confinement" and up to 10 5 by additionally 
shifting the third neighboring holes outwards [28]. The same approach has 
also been used to increase the Q-factor of a simple Ml-cavity [26J in diamond 
with one missing hole in the center from » 200 to 3 • 10 4 [29] by optimizing 
the next neighboring holes in the x- and y-direction and up to 7 • 10 4 [18] 




(1) 




[18j. 
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by fine-tuning the holes farther outwards. Together with a modal volume 
of V e g = 1.1 (X/n) 3 a theoretical Purcell-Factor of 4580 can be achieved by 
coupling a color center in diamond to the optimized cavity structure. 

In order to increase the Purcell-Factor even further, here we consider a 
cavity design with very small modal volume: By shifting two adjacent holes 
outwards a so-called MO-cavity is introduced [26J. Extending the design ap- 
proach of "gentle confinement" to remote holes and by comparing simulated 
and ideal field distributions in Fourier and real space, we show that the cavity 
Q-factor is improved significantly by optimizing the surrounding holes not 
only in the close vicinity of the point defect but also farther away. Exem- 
plarily, we consider in detail the design process of a MO-cavity in diamond 
and in silicon. 

Due to significant progress in diamond processing techniques, first exper- 
imental demonstrations of cavity modes in photonic crystal defect cavities in 
diamond have been achieved (25]. However, the measured quality factors in 
diamond are more than one order of magnitude smaller than the theoretically 
predicted values. This limitation might be due to scattering losses from the 
grain boundaries and subsequent emission from a nonsmooth surface. An- 
other loss mechanism is material absorption of the nano-crystalline diamond 
slab. In our simulations, we investigate the influence of material absorption 
on the quality factor. We show that this influence can be approximated by 
two simple models, replacing time-consuming simulations. 

One possibility to fabricate thin films with enhanced optical quality is 
the use of single- crystal diamond membranes. Recently, single-crystal free- 
standing membranes [3QM32] and waveguide structures [331 El] in diamond 
have been produced using focused ion beam milling. The use of such single- 
crystal diamond membranes with small absorption losses might pave the way 
for realization of high-Q photonic crystal microcavities. 

2. Cavity design 

Our work is focused on two dimensional photonic crystals consisting of a 
triangular lattice of air holes with a lattice constant a in a thin membrane 
of diamond suspended in air. The periodic structure gives rise to a photonic 
band gap for TE-like modes [HJ [22]. By introducing a defect, light can 
be localized in three dimensions within small volumes: in the horizontal 
plane, light is localized due to distributed Bragg reflection, and in the vertical 
direction due to total internal reflection. 
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2.1. Computation method 

The calculation of the near field patterns of the cavity mode was per- 
formed with a finite difference time domain algorithm (FDTD) [35], using a 
freely available software package with subpixel smoothing for increased accu- 
racy [36]. The simulated PC structure is a 27a x 27 a super cell with height 
8 a, surrounded by perfectly absorbing boundary conditions (PML) [3?]. To 
extract the resonance frequency as well as the cavity Q-factor, we use a filter 
diagonalization method [3E]- The mode volume V e s is calculated by 



^ = fe(r)\E(r)Jd*r 
maxe(r)|£'(r)| 2 

The integral of the electric energy is taken over the whole computational cell. 
To reduce the simulation time and the amount of stored data we implement 
mirror boundary conditions in the x, y and z directions. To assure that the 
employment of mirror boundary conditions does not effect the calculated Q- 
factors, we reanalyze the structure at the beginning and at the end of our 
design process by applying an even mirror symmetry to the z = plane only 
in order to select TE-like modes. Calculations are initially performed with 
a resolution of 32 points per lattice constant. In order to check for errors 
due to numerical discretization we also use a resolution of 40 and 50 points 



per lattice constant for all optimized structures of Sec. |2.3| but find similar 
Q-factors. 

The Q-factor calculated by the FDTD algorithm might vary with the 
chosen resolution of the simulation. Therefore, it is interesting to compare 
the FDTD results with an alternative method which calculates the Q-factor 
via the stored energy in the cavity and the radiated power [351 ESHS] . This 
second method permits to draw conclusions about further improvement of the 
cavity geometry by analyzing the wave vector components inside the "light 



cone" (see Sec. 2.3). If the defect cavity is surrounded by a sufficiently large 
number of air holes, the cavity lifetime is mainly limited by radiation losses 
in the vertical direction [39]. The radiated power in the far field is then 
determined by the 2D Fourier transforms FT2 of the near field components 
E x , E y , H x , H y in a plane S at a distance Az above the photonic crystal slab 
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[33 SI]: 



V 



J Idk x dk y (3) 



8A 2 P 

|fen|<fc 

with I = \FT 2 (H y ) + V FT 2 (E X )\ 2 + \FT 2 (H X ) - V FT 2 (E y )\ 2 , 

where r] is ^/io/co and / is the radiated intensity. The integral runs over all 
k- vectors inside the "light cone" \kn \ < k. The cavity near field is computed 
using FDTD simulations after 400 time-steps. After that time transient 
effects, resulting from a cutoff of our excitation pulse, have died away. The 
2D Fourier transforms of the near field pattern are taken at a plane S in 
a distance Az = A/2 j4T] above the photonic crystal slab. To compare the 
radiation losses of different cavity designs, we normalize the radiated intensity 
/ to the stored energy U in the cavity where U is given by [3D] : 



U = \j (e e(r)\E(r)\ 2 + » \H(r)\ 2 )d 3 r (4) 

The quality factor Q is defined by the stored energy U in the cavity divided 
by the total radiated Power P per cycle [4"D] : 



Q = »ji, (5) 

where u denotes the angular frequency of the cavity mode. To calculate the 
Q-factor, the integral Q of the energy stored between two planes S above and 
below the photonic slab as well as the integral in equation Q are performed 
as a discrete sum. Finally we take the time average of the radiated power P 
and the stored energy U over one period. This final time- averaging can be 
avoided by calculating the imaginary part as well as the real part of the fields 
instead of the real part only. It is sufficient to evaluate the complex fields at 
one instant in time, since the energy of complex fields does not change over 
one period. 

2.2. Starting point: simple MO- cavity 

As a basic cavity design, we choose a point defect structure with a very 
small modal volume: By shifting two adjacent holes along the x-direction 



by a distance d, a so called MO-cavity is created [26J . Figure |1 (a) shows the 
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Figure 1: Design and field components of the MO-cavity: (a) Nomenclature of the holes 
surrounding the defect, (b) E y - and (c) iJ z -component of the MO^-cavity mode. 



nomenclature of the MO-cavity design, which will be used in the following. 
As a first example, we consider a MO-cavity in a lossless diamond membrane 
with a refractive index of n — 2.4 [12]. The cavity design supports only 



one resonant mode, whose E y and i^-components are shown in Fig. 1(b) 



1(c), respectively. The quality factor of this simple MO-cavity depends on 
three parameters: the background radius R of air holes, the slab thickness 
h, and the shift d of the two holes creating the defect. By successively 
modifying these parameters in steps of 0.01a within ranges R e [0.25a, 0.29a], 
h G [0.9a, 0.95a] and d e [0.13a, 0.17a], we find Parameter Set M0 A , listed 
in Tab. [TJ yielding the highest quality factor of Q ~ 24000 of all simple M0- 
cavities with a modal volume V e s = 0.390 (A/n) 3 at a resonance frequency 
u = 0.3652 (27rc/a). The quality factor seems to be independent of the slab 
thickness h in this parameter range. Therefore, we keep h = 0.91a for the 
remainder of the paper. 



R(a) 


h(a) 


d(a) 


Q 


u(2irc/a) 




0.26 


0.91 


0.15 


24000 


0.3652 


0.390 



Table 1: Starting point for the design optimization of the MO-cavity (Parameter Set M0 A ) 



2. 3. Optimization of quality factor 

The Q- factor of the simple M0 A cavity is more than a factor of four smaller 
than equivalent designs in silicon membranes [26] due to the lower refractive 
index of diamond. In order to achieve comparable Q-factors in diamond, 
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the cavity design has to be optimized, without delocalizing the cavity mode, 
i.e. without increasing the modal volume. In two dimensional photonic 
crystals light can be localized in three dimensions due to Bragg reflection 
in the horizontal plane and by total internal reflection (TIR) in the vertical 
direction by the air cladding. However, only plane wave components with 
in-plane wave vectors \k»\ > k are guided within the slab by TIR whereas 
all modes with wave vectors inside the "light cone" \k\\\ < ko are radiated 
into the air cladding according to Snell's law [21] . These vertical radiation 
losses are crucial for photonic crystals in diamond because of the relatively 
low refractive index. The lower refractive index results in a higher mid-gap 
frequency of the band gap compared to silicon-based photonic crystals and 
therefore results in a larger corresponding light cone [22] ■ A suitable approach 
to minimize these radiation losses in the vertical direction is the method of 
"gentle confinement" [271 [28] : Fourier- and real-space analysis reveals that 
radiation losses can be reduced by tailoring the cavity mode profile to resume 
a Gaussian envelope function. In the case of photonic crystals, such "gentle 
confinement" can be obtained by adjusting the neighboring holes around the 
introduced defect, e.g. by reducing the radii or shifting some holes outwards. 
The design process of the MO-cavity is structured as follows: 

Step 1 Optimization for background radius R = 0.26a 

(a) As a first optimization step, we fine-tune the holes c, e, k and 6, g, 
f in the vicinity of the defect highlighted by the shaded region in 



Fig. 1(a) , for a background radius R = 0.26a, such that the cavity 
mode profile along the y- and x-axis, respectively, fits a Gaussian 
envelope. 

(b) Thereafter, we adjust the holes m, p and s, v, t farther out along 
the y- and x-axis, respectively. 

Step 2 Optimization for background radii R G [0.27a, 0.29a] 

(a) As a next step, we check whether the background radius R and 
thus the in-plane localization of the modified mode is still optimal. 
Therefore, we repeat step 1(a) for changed background radii R G 
[0.27a, 0.29a]. 

(b) Finally, we repeat step 1(b) for a background radius R = 0.28a. 
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(a) MO" 4 cavity mode 
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Figure 2: Comparison between the radiated intensity / normalized to the stored energy U 
of (a) a simple M0 A mode (b) a M0 B mode, where the holes along the y-axis are optimized 
and (c) a MO^ cavity mode where additionally the holes along the ir-axis are fine-tuned. 



Step l.a: Optimization of the next-neighboring holes for background radius 
R = 0.26a 

The modest quality factor Q of a simple M0 A cavity is explained by con- 
sidering the radiation losses in the vertical direction. The radiation intensity 
I of the M0 A cavity calculated according to equation ^ is shown in Fig. 
2(a) In order to properly compare radiation losses of different cavity de- 



signs, we normalize all plots of radiation intensity I to the energy U stored 
in the cavity. The white circle denotes the light cone boundary \k»\ < k . 
Note the logarithmic color scale to visualize weak field components. The 



red area inside the light cone in Fig. 2(a) indicates large leaky components: 



The radiated power inside the light cone divided by the stored energy is 
P/U ~ 7.9 ■ 10~ 5 (c/a). We start our design process by optimizing the cavity 
structure along the y-axis. In order to identify which holes have to be ad- 
justed, we analyze the field distribution along the y-axis and compare it to 



an ideal Gaussian envelope function. Figure 3(a) shows the ^-distribution 
of the M0 A cavity mode (black rectangles). The ideal mode profile, shown 
in red, is given by a sinusoidal fundamental wave multiplied by a Gaussian 



envelope function (green curve in Fig. 3(a)) |27J. The shaded regions indi- 
cate deviations of the cavity mode from the ideal profile. These deviations 
already start in the immediate vicinity of the cavity center, thus indicating 



the need for optimization of the holes c, e, k (see Fig. 1(a) ). 



As a first step, we vary the radii R c , R e and Rk of the holes c, e and k such 
that the cavity mode in the vicinity of the y-axis fits a Gaussian envelope. 
The optimal choice R c = 0.23a, R e = 0.24a, Rk = 0.25a yields a quality 
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Figure 3: Field amplitude of the MO" 4 - and M0 B -cavity mode along the y and x-axis, 
respectively: (a) i^-field of the MO^-cavity along the y-axis: The gray regions at the 
hole positions c, e, k indicate deviations of the calculated cavity mode (black) from the 
ideal mode profile (red) with a Gaussian envelope function (green) . (b) ^-distribution of 
the M0 B -cavity along the x-axis: The deviations from the ideal profile are important at 
positions b, s and t, whereas the deviations at position g are small. 



factor Q = 71400 with a modal volume of V e s = 0.390(A/n) 3 (Parameter Set 
M0 B in Tab. [2]). The 2D Fourier transform of the M0 B -cavity is shown in Fig. 



2(b) By adjusting the holes in the immediate vicinity of the defect the leaky 
components inside the light cone have been reduced to P/U ~ 2.8- 10 _5 (c/a). 



R(a) 


d(a) 


R c (a) 


-R e (a) 


R k (a) 


Q 


u(2irc/a) 


Kfr(A/n) 3 


0.26 


0.15 


0.23 


0.24 


0.25 


71400 


0.3622 


0.390 



Table 2: Optimization of first, second and third next-neighbor holes along the y-axis 
(Parameter Set M0 S ). 

As a next step, we repeat the procedure above for the x-direction by 
analyzing the ^-distribution along the x-axis as shown in Fig. 



3(b) The 



deviations from the ideal field distribution are important at the hole positions 
b, s and t, whereas the deviations at position g are quite small. The mismatch 
at the position of the next-neighbor holes b indicates, that the displacement 
d of the holes along the x-axis is no longer optimal, after we have optimized 
the holes in the vicinity of the ^/-direction. 

By slightly increasing the shift d = 0.16a, keeping all other parameters 
of the M0 B structure fixed, the quality factor can be improved significantly 
up to Q = 121500 with a modal volume V e g = 0.388(A/n) 3 . An additional 
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modification of the radius Rb of the next-neighbors does not lead to a further 



improvement of the cavity lifetime. As indicated in Fig. 3(b) the field distri- 
bution at position g seems to be already optimal. Indeed neither changing 
the radius R g nor shifting the holes g along the x-axis leads to an increase of 
the quality factor. Almost the same holds for the holes / farther away from 
the defect. By reducing the radius Rf the quality factor can be improved 
up to Q = 132000. The cavity design optimized in the x- and y-direction 
for background radius R = 0.26a is summarized in Parameter Set M0 C in 
Tab. [3] The radiation intensity / is shown in Fig. 2(c)[ The increase of the 
Q-factor compared to the M0 B structure is attributed to a reduction of the 
radiated power in the vertical direction to P/U ~ 1.6 • 10~ 5 (c/a). 
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a] 


d[a] 


R c [a] 


Re 


a] 


Rk 


a] 


R f 


a] 


Q 


u[2ttc/o\ 


Kff[(A/n) 3 ] 


0.26 


0.16 


0.23 


0.24 


0.25 


0.25 


132000 


0.3606 


0.388 



Table 3: Optimization of the holes along the x- and y-axis for background radius R = 0.26a 
(Parameter Set M0 C ). The radii not listed here correspond to the background radius R. 



Step l.b: Optimization of remote holes for R = 0.26a 

After the fine-tuning of the holes in the vicinity of the defect (see shaded 
region in Fig. 1(a) ), we consider the holes m, p and s, v, t at larger distances 
along the y- and respectively. Judging from the radiation intensity I 

of the M0 C cavity (Fig. 2(c) ), the wave vector components along the fc^-axis 



are already very small. They represent about 10% of the components in the 
light cone. Indeed, neither changing the radii R m nor R p of the holes m, p 
along the y-axis yields further improvement of the cavity Q. 

As a second step, we consider the holes s, v, t along the x-axis. In Fig. 
3(b)| deviations from the Gaussian envelope are visible at positions s and t. 
By slightly reducing the radii R s = 0.25a and R t = 0.25a the quality factor 
can be improved up to Q = 146500 with a modal volume V e s = 0.390(A/n) 3 . 
Additional modification of the holes v does not lead to further improvement 
of the cavity Q. The optimal choice of the hole radii R m , R p and R s , R v , 
R t in the y- and x-direction, respectively, are summarized in Parameter Set 
M0 D in Tab. |4| The other parameters not listed in Tab. [4] correspond to the 
parameters of the M0 -cavity. 

Step 2. a: Optimization of the next-neighbors for different background radii R 
At the beginning of our design process we adjusted the background ra- 
dius R for a simple MO-cavity. We now want to verify whether the choice 
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Rm 


a 




R p 


a 




R s 


a 




R v 


a 




Rt 


a 




Q 


u[27ic/a 




Kfr[(A/n) 3 ] 


0.26 


0.26 


0.25 


0.26 


0.25 


146500 


0.3604 


0.390 



Table 4: Optimization of the holes at larger distances from the defect for background 
radius R — 0.26a (Parameter Set M0 D ). The other parameters are the same as in Set 
M0 C . 



of R = 0.26a is still optimal. On the one hand, as discussed in Ref. [IB], 
larger background radii R lead to an increase of the width of the band gap 
and hence to a better in-plane mode confinement. On the other hand, an 
increase of the radius R also leads to larger wave vector components inside 
the light cone and thus to larger radiation losses in the vertical direction. 
Therefore, the choice of R is a tradeoff between these two opposed contri- 
butions. By optimizing the holes around the defect, vertical radiation losses 
have been reduced significantly. We now want to check whether the in-plane 
confinement of the mode can be improved as well, by increasing the back- 
ground radius R G [0.27a, 0.29a] in steps of 0.01a. Using the very same 
design process as described in step l.a, we fine-tuned the holes c, e, k and 
6, g, f in the y- and x-direction, respectively, for every background radius 
R G [0.27a, 0.29a]. The highest quality factor of an optimized cavity geom- 
etry is obtained for background radius R = 0.28a (Parameter Set M0 E in 
Tab. [5]). For an optimized cavity structure, the confinement of the mode 
in the in-plane direction is considerably improved for a background radius 
R = 0.28a leading to larger quality factor Q = 226600 and a smaller mode 
volume Ves = 0.350(A/n) 3 . 



R[a] 


d[a] 


R c 


a] 


Re 


a] 


Rk 


a) 


R f 


a] 


Q 


u[27ic/a] 


v cS [^M 3 ] 


0.28 


0.16 


0.22 


0.24 


0.26 


0.26 


226600 


0.3673 


0.350 



Table 5: Optimized MO-cavity for background radii R = 0.28a (Parameter Set M0 B ). 



Step 2.b: Optimization of remote holes for background radius R = 0.28a 

After the optimization of the holes in the vicinity of the defect (see shaded 
region in Fig. 1(a) ), we again investigate the influence of the remote holes m, 
p and s, v, t. Fig. 4(a) shows the normalized radiated intensity of the MO^ 
cavity mode: Wave vector components are still left inside the light cone. The 
radiated Power P normalized to the stored energy is P/U ~ 1 • 10 _5 (c/a). 
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Figure 4: Comparison between the radiated intensity normalized to the stored energy for 
background radius R = 0.28a of (a) a M0 E mode (b) a M0 F mode, where the holes both 
along the x- and y-axis are optimized, (c) E y -field of the M0 F -cavity mode along x-axis: 
The mode profile fits well to a Gaussian envelope. 



This suggests, that further improvement of the high-Q MO^-cavity can be 
obtained by adjusting the holes at larger distances from the defect. 

Like before, we start to optimize the holes m, p along the y-axis, which 
leads to only slight improvement of the Q-factor to 228500 for radii R m = 
0.31a, R p = 0.34a. As a final step, we consider the holes s, v, t along the x- 
axis (see Fig. 1(a)). The optimal choice of the radii R s = 0.29a, R v = 0.27a, 
R t = 0.28a yields another significant decrease of the leaky components inside 

The reduction of radiation losses to 



4(b) 



the light cone, shown in Fig. 
P/U ~ 6.9 • 10~ 6 (c/a) leads to a considerable increase of the quality factor 
up to Q = 320000 with a modal volume V e s = 0.350(A/n) 3 (see Parameter 
Set M0 F in Tab. [6]). Coupling a color center to such an optimized cavity 
leads to a theoretical Purcell-Factor of F = 69500. The ^-component of the 
M0 F -cavity along the x-axis is shown in Fig. 4(c) After carefully fine-tuning 



the holes around the defect the mode profile of the optimized structure fits 
well to a Gaussian envelope. 
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R P 
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Rt 
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V eS [(\/nf\ 




0.31 


0.34 


0.29 


0.27 


0.28 


320000 


0.3672 


0.350 



Table 6: Optimization of the holes along the x-axis (Parameter Set MO^). The other 
parameters are the same as in Set M0 E . 
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Figure 5: Fabrication tolerance test: Dependence of (a) the Q-factor and (b) the modal 
volume VcB of the optimized MQ F cavity on the radius R c and the displacement d. 



2.4- Fabrication tolerance tests 

In view of experimental realization of photonic crystals in diamond, it is 
interesting to investigate the influence of the cavity Q on fabrication toler- 
ances. As discussed in section 2^3 the quality factor is significantly enhanced 
up to Q = 320000 by modification of the holes surrounding the defect. In 
this section, we investigate how sensitive the high-Q modes are to variations 
of the neighboring holes. Exemplarily, the next-neighboring holes b and c are 
considered. Starting from the optimized M0 F -cavity (see Tab. |6f, we vary 
the displacement d G [0.15a, 0.17a] and the radius R c G [0.21a, 0.23a] in 
steps of 0.01a. The dependence of the quality factor and the modal volume 
on the surrounding structural parameters are shown in Fig. [5| The displace- 
ment d seems to be a crucial parameter. If d changes from 0.16a to 0.15a, the 
Q factor drops down to 4 — 6 ■ 10 4 , for all radii R c . In contrast, if we keep the 
shift d = 0.16a of the holes b fixed and vary the radius R c G [0.21a, 0.23a] 
only, the Q-factor of the M0 F cavity always exceeds 10 5 . This is favorable 
considering that it is easier to fix the central position of the air holes than 
to determine their radius using common etching techniques. If we envisage 
a central wavelength A = 637 nm of the M0 F -cavity equal to the emission 
wavelength of the NV-center, the lattice constant is a = 234 nm. To achieve 
large quality factors in practice, high precision structuring techniques with 
fabrication tolerances smaller than 3 nm are required. One available tech- 
nique is to use a focused beam of gallium ions (energy 30 keV) to pattern the 
diamond membrane. The spot size of the gallium ion beam is in the order of 
several nanometers, so structures can be fabricated with a resolution of 6 nm 
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Another technique which has already been used to fabricate photonic 
crystals in diamond (25] is reactive ion etching achieving radii of 80 nm ± 
5 nm jH] . The radii of the milled air holes are sensitive to the dose used in 
the e-beam lithography to pattern the etch mask and to the applied etching 
parameters. 

In contrast to the cavity Q, the modal volume seems to be nearly inde- 
pendent of the fine-tuning of the surrounding holes. As shown in Fig. 5(b)| 
the modal volume varies between V e s = (0.336 — 0.358) (A/n) 3 . This is a huge 
advantage of the MO cavity design in view of achieving large enhancement of 
the spontaneous emission. Even if the quality factor is reduced to Q ~ 4 • 10 4 
due to fabrication tolerances, theoretical Purcell factors of F ~ 8500 can be 
obtained by coupling an emitter to the M0 F -cavity. 



2.5. M0 cavity design in a silicon slab 

As discussed above, we have optimized the MO-cavity structure in dia- 
mond such that vertical radiation losses have been reduced significantly yield- 
ing a Q-factor even larger than predicted for MO-cavity structures in silicon 
|26j . Because of the larger refractive index of silicon, this suggest that further 
improvement of the MO-cavity geometry can also be obtained in a silicon slab 
using the same design process as described above. As a starting point for 
the optimization in silicon we choose the MO-cavity structure presented in 
Ref. [26] yielding a quality factor of Q ~ 135000 with a modal volume of 
V e s = 0.288(A/-n,) 3 . After carefully fine-tuning the background radius R and 
the thickness h of the slab as well as the surrounding holes according to the 
procedure described in section 2.3, we find the optimal M0 -structure listed 



in Tab. [7j The other radii not listed in Tab. [7] correspond to the background 
radius R = 0.26a. As a result, we gain an improvement in cavity Q by more 
then a factor of three with Q « 458000 and V cS = 0.260(A/n) 3 . This result 
confirms that adaption of air holes even far away from the defect generally 
improves the cavity lifetime significantly. 
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0.26 


0.8 


0.14 


0.24 


0.24 


0.25 


0.24 


0.28 


0.25 


458000 



Table 7: Optimization of the MO-cavity in a silicon slab (Parameter Set M0 Sl ). The radii 
of the air holes not listed here are equal to the background radius R. 
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Figure 6: Comparison between the Q-factors calculated from the radiated power according 
to equation ^ and the FDTD simulations. On the one hand the radiated power is 
determined by the Poynting vector (red points) and on the other hand by integrating over 
all components inside the light cone using equation Q (green triangles). 



3. Radiation Q-factors versus FDTD results 



In section |2.3 we calculate the radiated intensity from the 2D spatial 



Fourier transforms of the cavity near field. The leaky components inside 
the light cone indicate in which direction further improvement of the cavity 
geometry is necessary. Furthermore, the radiated intensity I gives not only 
qualitative information about the radiation losses but also permits to draw 
quantitative conclusions. By summing over all wave vector components in- 
side the light cone according to equation Q, we can calculate the radiated 
power P. This offers in alternative method to determine the cavity Q-factor 
according to equation (|5]): dividing the stored energy U in the cavity by the 
vertical radiated power P and multiplying by the resonance frequency u. In 
this section, we compare the radiation Q-factors with the FDTD results. The 
Q-factors calculated by the FDTD algorithm might depend on the chosen res- 
olution (black rectangles in Fig. [6]). The errorbars in Fig. [6] (exemplarily for 
structures M0 A , M0 E , M0 F ) show the variation of the FDTD results with 
the change of resolution from 32 to 50 points per lattice constant. Moreover, 
the application of mirror symmetries to reduce the computational cell, might 
also cause an overestimation of the Q-factors [5H]. Therefore, it is interesting 
to have a simple alternative method to determine the Q-factors in order to 
identify numerical artifacts. 

We calculate the near field patterns as well as the stored energy U in the 
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cavity using FDTD simulations. The vertical radiation losses are determined 



via the Fourier transforms of the near field (see Sec. 2.1 ). The spatial Fourier 
transforms are calculated in a plane S at position z — A/2. To determine the 
vertical radiated power, the integral Q over all wave vector components in- 
side the light cone is performed as a discrete sum. Fig. [6] shows a comparison 
between the Q-factors calculated via radiated power (green triangles) and the 
FDTD results (black rectangles) for the cavity designs in diamond presented 



in sections 2J2 and 2J3 This simple calculation procedure reproduces well 
the FDTD results for different cavity designs. For structures M0 A , M0 B , 
M0 C and M0 D with small background radius R = 0.26a, deviations from the 
FDTD results are maximally 20% (structure M0 D ). These deviations are 
mainly due to the fact, that the resonance frequency is close to the upper 
edge of the photonic band gap leading to imperfect lateral confinement of 
the cavity mode. Albeit, the radiation Q-factor, taking into account verti- 
cal radiation losses only, reflects well the improvement of the cavity lifetime 
by optimizing the photonic crystal geometry. For structures M0 E and M0 F 
with enlarged background radii R = 0.28a, the deviations of the radiation 
Q-factor from the FDTD results are smaller than 6%. With increasing back- 
ground radii, the photonic band gap increases and the resonance frequency 
of the cavity mode lies deeper inside the photonic band gap resulting in a 
better lateral mode confinement. Therefore, by calculating the vertical ra- 
diation losses via Fourier transform of the near field, the radiation Q-factor 
of photonic crystal geometries with enlarged background radii is a good ap- 
proximation to the FDTD results. 

A third approach for determining the Q-factor uses the radiated power 
as well as the stored energy obtained from FDTD simulations. The total 
radiated power in this case is computed by the integral of the Poynting 

—* — * 

vector y\z(E* x H) in all three directions through planes at the edge of the 
simulation cell. The associated Q-f actors for cavity geometries in diamond 
are shown by the red points in Fig. [6] The agreement with the FDTD results 
is very well. 

Calculating the stored energy as well as the radiated power, either using 
FDTD simulations or via spatial Fourier transforms of the near field, offers 
an alternative way to determine the cavity Q-factor. This method can be 
used to check the results of the FDTD simulations, usually obtained from 
the analysis of the cavity temporal decay, and to easily identify possible 
numerical artifacts. 
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4. Absorption 

For the practical realization of photonic crystal microcavities in diamond 
it is essential to investigate the influence of the material absorption on the 
cavity Q-factor. To achieve effective optical coupling in practice, ultrapure, 
high-quality diamond films are required. Actually, our measurements show 
|45j . that nano-crystalline diamond films strongly absorb light at the emis- 
sion wavelengths of the NV, SiV and NE8 color center (~ 630 nm - 800 nm) 
due to graphite or amorphous carbon in the grain boundaries [46] . The ab- 
sorption coefficient is determined by measuring the transmission through a 
nano-crystalline diamond membrane (p-BeSt, Innsbruck, Austria). The re- 
duction of the transmission signal due to surface roughness (15 nm rms) is 
estimated to 2%, while the material absorption leads to a reduction of the 
signal of 11% (for a film of 300 nm thickness). Therefore, we assume that ma- 
terial absorption of the dielectric background is the primary loss mechanism 
of microcavities in photonic crystals in nano-crystalline diamond films. 

To investigate the influence of material absorption on the cavity Q-factor, 
we implement artificial resonances in the dielectric function such that the 
imaginary part 6j yields the desired resonant absorption whereas the real 
part e r remains almost constant [IB]. The results of the FDTD simulations 
(black rectangles in Fig. [7]) reveal that the Q-factor strongly depends on 
the absorption coefficient a of the diamond material. For our ultra-nano- 
crystalline diamond films we measure maximum absorption coefficients a ~ 
4000 cm" 1 at the emission wavelength A = 637 nm of the NV center |45j. 
Starting from a theoretical quality factor Q 50000 of a MO-cavity with 
optimized next neighboring holes (R = 0.26a, d = 0.15a, R c = R e = 0.24a, 
h = 0.91a), the material absorption of the diamond film would reduce the 
quality factor to only Q < 100. Nevertheless, coupling an emitter to a 
MO-cavity, a Purcell-Factor of about F = 15 is still achievable because of 
the extraordinary small modal volume. Due to recent progress in growing 
high-quality diamond films by chemical- vapor deposition [17] , absorption co- 
efficients of a ~ 200 cm -1 [IB] can be obtained. With such low loss diamond 
membranes an experimental cavity Q-factor of 1150 can be achieved, yielding 
a Purcell-Factor of F = 220. 

Calculation of the quality factors for different absorption coefficients a 
using FDTD simulations is a very time-consuming task. Therefore, an ana- 
lytical approximation to evaluate the Q-factor for lossy materials would be 
very helpful. The dependence of Q on a can be described by a simple model: 
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Absorption coefficient a (cm 1 ) 

Figure 7: Influence of the absorption coefficient a of the diamond material on the Q-factor 
of the MO-cavity (R = 0.26a, d = 0.15a, R c = R e = 0.24a, h = 0.91a): The FDTD results 
are shown in black. The dependance can be described by a simple model of a linear cavity 
filled with a homogenous diamond slab or filled with the actual photonic crystal structure. 

In the case of a linear cavity with a homogenous diamond slab between two 
mirrors the total quality factor Q is given by [191 150] : 

Q- 1 = QcT 1 + q:L (6) 

where Qq denotes the quality factor of an ideal, lossless cavity and Qabs the 
quality factor due to material absorption. As described in Refs. [18j 150] . 
for weakly absorbing material e r 3> the absorption quality factor can be 
written as Q a b s = k/a, where the absorption coefficient is given by a = 
koei/y/e^. The wave vector in the dielectric is k — y/e^ko, where ko denotes 
the wave vector in vacuum. Under the assumption that most of the electric 
field is concentrated in the dielectric, the material losses are given by: 

Qate = k/a = e r /ei (7) 

The blue curve in Fig. [7] shows that this first approximation predicts well the 
dependance of the cavity Q-factor. The deviations from the FDTD results 
are smaller than 13%. Actually, equation ([7]) denotes a lower limit to the 
FDTD results for lossy materials, because the electric field is not completely 
concentrated within the dielectric but also leaks into the air holes. 

In a more general approach, we want to take into account the actual pho- 
tonic crystal structure. For weakly absorbing material e r 3> the imaginary 
part 6j can be considered a small perturbation Ae = zej of the dielectric 
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function e(f). This perturbation Ae results in a first-order correction of the 
resonance frequency Au = iuji (5TJ [52] : 



A „ = uJ r jAe(f)\E(r)\ 2 d 3 r 
2 j er (f) \E(r)\ 2 d 3 r ' 

where u r denotes the resonance frequency of the ideal, lossless cavity struc- 
ture. Considering that the Q-factor is defined as Q = — the absorption 
quality factor can be written as: 

_ j e r {r)\E{r)\ 2 d 3 r 
Je t (f) \E(r)\ 2 d 3 r 

The integral in equation ^ is taken over the entire photonic crystal slab 
with the electric field E{r) of the ideal, lossless photonic crystal structure 
as a weighting factor. Therefore the electric field has to be calculated only 
once for an ideal cavity using numerical methods such as FDTD simulations. 
Afterwards it can be reused to calculate Q' abs for different lossy materials. The 
overlap integral in equation ^ of the electric field and the dielectric function 
is performed as a discrete sum. The red curve in Fig. [7] shows the absorption 
quality factor Q'^ calculated for different absorption coefficients a of the 
dielectric background. The deviations from the FDTD results are mostly 
due to different procedures used to calculate the overlap integral in equation 
The discrete summation leads to a slightly higher Q- factor of about 
11% compared to an exact integration. We estimate the error by calculating 
the effective dielectric constant J e(r)\E(r)\ 2 d 3 r / J \E(r)\ 2 d 3 r performing the 
integrals as a discrete sum and comparing it to FDTD simulations. 

Both simple models obtained above predict well the dependence of the 
quality factor for lossy materials. These simple equations permit to calculate 
an upper and a lower limit to the FDTD results without relying on time- 
consuming simulations. 



5. Conclusion 

We have presented an optimized MO-microcavity design in a diamond- 
based photonic crystal. The small defect is introduced by shifting two 
holes outwards along the x-axis yielding an extremely small modal volume 
of V e ff = 0.35(A/n) 3 . Using Fourier- and real-space analysis, we improved 
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the cavity structure by systematically varying the radii of the neighboring 
holes around the defect. The improved design yields a quality factor of 
Q = 320000. Coupling a color center to such an optimized cavity structure 
leads to a theoretical Purcell-Factor of F — 69500, which is one of the largest 
enhancement of the spontaneous emission that has been predicted so far for 
photonic crystal point defects in diamond. Using the same design process, we 
have additionally optimized the MO-cavity in silicon gaining a factor of three 
in cavity Q. The significant improvement of the cavity lifetime by "gentle 
confinement" shows, that it is worth looking closer at the surrounding holes 
even far away from the defect. 

In view of the practical realization of photonic crystals in nano-crystalline 
diamond films, material absorption of the dielectric background has to be 
taken into account. The influence of absorption losses on the cavity Q-factor 
can be described by a simple model of a linear cavity: In first approximation, 
we consider a homogenous diamond slab between two perfect mirrors. This 
simple description denotes a lower limit to the FDTD results, because the 
electric field is not completely concentrated within the dielectric, but also 
leaks into the air holes. The second approximation takes into account the 
actual photonic crystal structure. These two simple models offer an efficient 
way to predict quality factors of photonic crystal cavities for lossy materials 
without relying on time-consuming FDTD simulations. Whether this simple 
model of a linear cavity is valid for a larger range of defect geometries in 
photonic crystals is subject of current work. Preliminary results show that 
equation §f§ predicts the absorption quality factor Q a bs even better for large 
photonic crystal point defects. 

One possibility to overcome the limitation of the Q-factor by material 
absorption is the use of single crystal diamond films. Recent progress in 
fabrication of thin free-standing single crystal diamond membranes, might 
pave the way for efficient direct coupling of color centers to high-Q diamond- 
based photonic crystal microcavities. 
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